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Abstract We have studied dimerized spin systems by 
realizing the cluster expansion to high order. We have 
extended our previous dinrer expansion for one-dimen- 
sional systems to cover weakly interacting chains for a 
quantitative description of three dimensional materials 
like (C4Hi2N2)Cu2Cl6 (= PHCC) and KCuClg. By com- 
parison with recent inelastic neutron scattering data we 
are able to determine the exchange energies between 
individual spins. We have further investigated the in- 
commensurate region of zigzag chains with isotropic ex- 
change coupling constants near the disorder-line where 
the dispersion curve exhibits a minimum at a finite wave- 
vector. Our approach clearly shows the gradual transi- 
tion between the minimum of the dispersion at wavevec- 
tor and wavevector n within this region. The extent 
of the incommensurate regime is given analytically in an 
expansion in the coupling constants. 
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1 Introduction 

Spin systems consisting of dimers as basic building blocks 
•■ describe a large number of real materials when additional 
' exchange couplings of smaller magnitude between dimers 
■ are considered. For a quantitative treatment of the ele- 
' mentary excitations of such systems we have developed 
the dimer perturbation theory [nlH| to cover dimers in- 
teracting in two and three spatial dimensions up to high 
orders in the interdinrer interactions. In this contribu- 
tion we present the application of this method to study 
quantitatively the elementary excitations in PHCC and 
KCuCla (sect, y) and to determine the incommensurate 
regime of the one-dimensional zig-zag chain close to the 
dimer point (sect. |). The resuhs for PHCC and KCuCla 
were motivated by and will be compared to spectra ob- 
tained in recent inelastic neutron scattering experiments. 
For weakly interacting dimers the cluster expansion 
is a powerful tool to calculate power series for properties 
like groundstate energy p] or one magnon excitations B . 



Due to the large number of different exchange couplings 
the obtained power series have one parameter only which 
is an overall scaling of the exchange couplings and finally 
set to one. Different sets of exchange parameters require 
new calculations whereas the dependence in q-space is 
given completely. The power series which are correct for 
the infinite system are build from finite clusters. In our 
case going to sixth order the number of topological dif- 
ferent colored clusters is 18084 for KCuCls and to fourth 
order 405 for PHCC. 

In our notation the basic intradimer interaction (which 
sets the energy scale) is denoted by </; interdimer inter- 
actions are denoted by Jn^„) with [Iran) giving the di- 
rection between the dimer centers (in units of the lattice 
constants) and [ij) specifying the pair of spins which in- 



teract with exchange energy JY,'" s. The z— component 

2 on a given 
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{Ivan 

1 to spin i 



of the vector from spin i 
dimer is taken to be > 0. 

Using this method, it is possible to determine the ex- 
change interactions between the individual spins fornring 
the dimers whereas the standard method to analyze sys- 
tems of interacting dimers is an RPA-like approach Q| 
which allows to deternrine only the effective dimer inter- 

0/2. 



actions J?- = (J^HL + J^^L - J^.. - J^J'' 



{Imn) V (Imn) {Irrtn) {Iran) '^ {Imn) 

This effective approach reproduces the leading terms of 
the full series; the full series, however, has additional 
terms starting in second order which turn out to be inr- 
portant. 



2 Chains interacting in two and three 
dimensions: PHCC and KCuCla 

In this section we present our results for the real interact- 
ing dimer conrpounds (C4Hi2N2)Cu2Cl6 (Piperazinium 
Hexachlorodicuprate — PHCC) and KCuCla. In both 
these materials Cu^+— ions with S = ^ are responsible 
for the magnetic properties. 

PHCC was investigated recently by inelastic neutron 
scattering J6| and found to be effectively two-dimensional 
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(i.e. J(i^j^) = for m ^ 0). The compound can be 
considered as built from interacting spin ladders. The 
inelastic neutron scattering spectra were analyzed us- 
ing the effective dimer model. Satisfactory agreement 
with the data was obtained when effective interactions 
Almn) ^^^^ [Imn) — (200) and (002) were introduced as 
nonzero; however, no obvious exchange paths exist for 
these separations. In the following we present an analy- 
sis of the spectra with interactions Jn^j^) between neigh- 
bouring dimers only. The results shown in figure are 
obtained with the set of exchange interactions given in 
table |l|. A detailed analysis shows that for PHCC the 
additional contributions to the dimer scries which are 
not taken into account by the effective dimer approach 
are essential to lower the dispersion curves. On the other 
hand we find from a comparison of the complete and the 
fourth order result of the effective approach that higher 
orders can be safely neglected. 
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Fig. 2 Dispersion of KCuCls for two directions: solid lines 
denote the result obtained from cluster expansion, circles are 
data from S . 
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Fig. 1 Dispersion of PHCC for two directions; solid lines 
denote the result obtained from cluster expansion, circles are 
data taken from H. 



The material KCuCla is similar in structure to PHCC. 
It contains, however, interdimer interactions also in the 
remaining (third) spatial direction. Inelastic neutron scat- 
tering experiments on this material have been published 
by two groups H,p| and have been interpreted using the 
effective dimer model. In earlier work [BJ we analyzed 
these data using dimer series expansions to fourth or- 
der. Recently we improved the computer implementation 
and in figure we show a spectrum covering all contri- 
butions up to sixth order. This result demonstrates the 
convergence of our series expansion for the exchange con- 
stants between the individual spins relevant for KCuCla 
as given in table |l|. It confirms our previous result that 
the one-dimensional subunits in this material are chain- 
like structures with alternating couplings rather than 
ladders. 



3 Incommensurate regime in zig-zag chains 

In this section we consider a one-din'^ensionaj array ^f 



dimers forming an interacting zig-zag chain, J 



J, 



(22) 
(100) 



J2, J, 



(21) 
(100) 



(11) 
(100) 



Ji. This system has the inter- 



esting feature of an incommensurate regime where the 
minimum of the basic magnon dispersion curve is at fi- 
nite wavevector q = qmin, < qmin < tt. The incom- 
mensurate regime is close to the disorder line (Shastry- 
Suthcrland line) Ji = 2 J2 , where noninteracting dimers 
form the exact ground state, but requires slightly larger 
values of J2 for given Ji. Using the dimer expansion to 
sixth order we have determined the limiting lines of the 
incommensurate regime as 

Jld(J2) <Jl< Jlu{J2) 

Jid{J2) = 2J2 - 2J2 + 2J| - ^ J^ + ^4 + 0{Jl) 



Jlu{J2) — 2J2 — 2J2 + 2J2 



2 ^ 



^J^ + OiJ^) (I) 



Thus the deviation from the Shastry-Sutherland line is 
of second order in the interdimer interaction whereas 
the width of the incommensurate regime is of fourth or- 
der. The incommensurate regime is shown in figure in 
a comparison of the results from high order expansions 
(which agree with the findings from the diagonalization 
of finite systems using the Lanczos algorithm) and the 
approximation of eqs. (|l|). The behaviour of the wave 
vector qmim corresponding to the minimum excitation 
frequency is also shown for one path through the incom- 
mensurate regime in twelfth order (thus improving on 
the results shown in ref. [^). 

Evidently, the incommensurate regime is tiny and 
therefore so far mostly of academic interest; however, if 
a material becomes available which allows tuning the ex- 
change parameters by e.g. alloying or applying external 
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pressure, our results provide a quantitative guide where 
to expect this regime of considerable qualitative interest. 
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Table 1 Considered interaction constants of KCuCls and 
PHCC. The energy scale is given by J{ooo)- 
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Fig. 3 Incommensurate regime: The solid line denotes the 
Shastry-Sutherland line, the other lines are eqs. (hi) up to fifth 
order. The inset shows the wavevector qmin for the minimum 
value of the dispersion on the line J\ — |(1.2 — J2). All 
symbols are obtained from expansions up to twelfth order. 



4 Conclusions 

Using the dimer expansion after computer implementa- 
tion we have investigated the low energy dynamics of 
the real materials PHCC and KCuCla and of the zig-zag 
chain, a model of theoretical interest in its incommensu- 
rate regime. The extent of the incommensurate regime 
in the zig-zag chain is satisfactorily obtained in an ana- 
lytical approximation close to the dimer point and this 
result provides simple quantitative predictions for pos- 
sible candidates for this interesting phase. The magnon 
spectra of PHCC and KCuCla as measured in inelas- 
tic neutron scattering experiments can be used to de- 
termine the exchange energies between the individual 
spins forming the basic dimers by applying the dimer 
expansion to fourth, resp. sixth order. KCuCls turns 
out to be appropriately described as a system of in- 
teracting chainlike structures with alternating couplings 
rather than ladders. For PHCC a consistent set of ex- 
change energies is presented for interdimer interactions 
between dimers which are nearest neighbours only. For 
the material TlCuCla, which is similar to, but stronger 
interacting than KCuCla, the application of the dimer 
expansion approach to the magnon spectra will be pre- 
sented together with new experimental results [lO||. 
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